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Abstract 

The quantum average energy decay and the purity decay are studied for a system 
particle as a function of the number of constituents of a discrete bath model. The 
system particle is subjected to two distinct physical situations: the harmonic os- 
cillator (HO) and the Morse potential. The environment (bath) is composed by a 
finite number N of uncoupled HOs, characterizing the structured bath, which in 
the limit — t- oo is assumed to have an ohmic, sub-ohmic or super-ohmic spectral 
density. For very low values of N the mean energy and purity remain constant in 
time but starts to decay for intermediate values (10 ^ N < 20), where two dis- 
tinct time regimes are observed: two exponential decays for relatively short times 
and a power-law decay for larger times. In this interval of N decoherence occurs 
for short times and a non-Markovian dynamics is expected for larger times. When 
N increases, energy and coherence decay very fast and a Markovian dynamics is 
expected to occur. Wave packet dynamics is used to determine the evolution of the 
particle inside the system potentials. 
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1 Introduction 



The quantum dissipation and decoherence are analyzed here for an open sys- 
tem interacting with its environment by colhsion processes. The problem 
(System+Environment+Interaction) is conservative but, due to energy ex- 
change between system and environment, the system can be seen as an open 
system with dissipation. Such theoretical model has been proposed and stud- 
ied in quantum [TP|5] and classical PI^IH] systems. Various models were de- 
veloped to treat such open system. We start by mention the methods which 
have focused on an explicit quantum dynamical treatment of the system-|-bath 
dynamics: the multiconfiguration time-dependent Hartree (MCTDH) tech- 
nique for wave packet propagation [5], the Gaussian-MCTDH approach [6], 
the effective-mode representation [7], and the local coherent-state approxima- 
tion to system-bath dynamics [S]. Another way of studying the system+bath 
dynamics is to solve the non-Markovian master equations [9|T0] . including 
some semiclassical approaches [TT|12|13] . 




Here a non-Markovian quantum trajectory theory is used, named non-Mar- 
kovian quantum-state diffusion, that describes the dynamics of a quantum 
system coupled linearly via position to a environment [T^15fl6|[T7] . Differ- 
ently from most the previous studies, the environment is composed by a fi- 
nite number, A^, of uncoupled HOs, as studied before for classical systems 
[THUn]. In such cases we say to have a discrete or structured bath. Structured 
baths have become important since they describe realistic physical situations 
of non-equilibrium physics. One major recent application is the energy transfer 
between a light-harvesting protein and a reaction center protein [20]. Under- 
standing the physics encountered in this process will help exploring its huge 
potential [2T]. Some experiments in which the non-Markovian behavior arises 
due to the discreteness of the bath can be mentioned: high-Q microwave cav- 
ities, quantum optics in materials with a photonic band gap, output coupling 
from a Bose-Einstein condensate to create an atom laser [22]|23f24f25ll26j or the 
decoherence phenomenon [27ll28f29|l30] . From the classical point of view, in the 
context of finite baths with uncoupled HO, there are some works which analyze 
the effect of discrete (structured) baths on the system energy decay [l|18fl9] . 
The main result found is that finite baths may induce a non-Markovian dy- 
namics on the system particle. From the quantum point of view, most works 



[Hl6|8l[T3|30|31] focus on the analysis of changing the frequencies distribution 

and the coupling strength between system and bath. 

We analyze systematically the effect of increasing the number of oscillators 
from the bath on the system energy decay. In most cases we vary A^ = 
1 — >■ 100. As A^ increases, the quantum simulations, usually done here for 500 
realizations of the bath, need very long computational times. For very low 
values of A^ the mean time energy decay keeps constant but starts to decay 
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for intermediate values of A^. For 10 < < 50 a transition to exponential 
decay for shorter times and a power-law decay for larger times is observed. 
For higher values > 100 only exponential decay is observed. In the interval 
20 < < 50 an additional power-law decay between the two previously 
mentioned behaviors is observed. 

The paper is organized as follows. In Section [2] we introduce the basic con- 
cepts of non-Markovian quantum-state diffusion and derive a new simplified 
stochastic Schrodinger equation valid for finite baths. In Section [3] the simpli- 
fied equation is applied to analyze the A^ dependence of some useful physical 
quantities such as the energy decay, the average position, and the purity (de- 
coherence). In this section harmonic potential is considered. In Section |4] a 
Morse potential is used for the system particle. In the last section we present 
a summary and the conclusions. 



2 Linear non-Markovian quantum- state diffusion 



The non-Markovian quantum-state diffusion (QSD) equation is based on a 
standard model of open system dynamics: a quantum system interacting with 
a bosonic environment with the total Hamiltonian 



Htot = H + ^gx{LaX + L+Oa) + hY^uJxatax, (1) 

A A 

where H is the Hamiltonian of the system of interest, L is a system operator 
coupling to environment and , cl\ are the raising and lowering operators. 



with the property [0^,^^] 
is given by 



The linear non-Markovian QSD equation 



-iH' iV^t) + Lzl \il)t) -L+ dsKit - s 



(2) 



where was assumed a factorized total initial state \ipo) |0i) IO2) ■ ■ ■ |0a) ■ ■ ■ , with 
an arbitrary system state \ijjo) and all environmental oscillators in their ground 
state |0a). 



zt = -^Y.9l4e'^'' (3) 

A 

is a complex Gaussian process with zero mean and correlations M [2:^*23] = 
K{t — s) = Y,x \g\\'^(^~^^^^^~^\ where M[] denotes the ensemble average over 
the classical driving noise. The form of K{t — s) corresponds to the zero- 
temperature limit. A is the number of oscillators from the bath, gx is the 
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coupling to oscillator A, ujx is its frequency and z\ = {xx + iyx)/\/2, where 
{xx,yx) are normal (Gaussian) distributed real number with zero mean and 
deviation one. The many realizations to solve the stochastic Schrodinger equa- 
tion are done over these normal distributed numbers. The functional deriva- 
tive of equation (j2]) can be written in the first order approximation as a linear 
function of \ipt) 



^ = 0{t, s, z*) \^|;,) = e-^(*-)Le^^(*-) . (4) 

Now, taking L = = q, which is the position operator from the system, the 
stochastic Schrodinger equation becomes 



h 



iPt) = -iH' l^i) + qzl l^i) - q dsK{t - s)e-^^(*-^)ge^^(*-^) . (5) 



In equation (|2]) and ([5]), the Hamiltonian H' = H{q,p) + q'^A{t), with A{t) = 

I 1 2 

J2x ^^^{cos uxt — 1), contains an additional potential term that turns out to 
be counterbalanced by a similar term arising from the memory integral [35] . 
Our approach is consistent with the Redfield theory [37] in the long-time limit 
when /q is replaced by /q°°. 

Using \tpt) = E„Cn(t) \(pn), H\(t)n) = e„|0n), Q = J2n,m(!nm\<Pn) < 0m|, and 

orthonormal condition < 0n|0n' >= ^nn', equation (|5]) transforms to 



hCn{t) = -ienCnit) - iA(t) ^ qnmqmm' Cm' (t) + Z* ^ qnmCm{t) 

mm' rn 

— X] 1nm'0mm'{t)Cm'{t), (6) 



where 



Q-i[i^X + {<^rn-<irn')]t — I 

Omm'{t) =< (Pm\0{t) \(j)m') = Qmm' Y] I^'aI^— 7 Tf, (7) 

^ -l[UJx + {em - £m'>\ 



and 



0{t) = f dsK{t - s)e-'^(*-^)ge'^(*-^). (8) 

JO 

For simplicity, we choose the Planck constant h = 1 throughout the paper. At 
next we analyze the following quantities: The quantum average of the system 
energy which is given by 
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(9) 



n 



and the quantum average of the system position 



{(l)t = (^tk l^i) = C*n{t)Cm{t)q, 



(10) 



and finally the purity which gives informations about the decoherence occurred 
in the system [T3] : 



where ps(t) = M [\ipt) {ipt\] is the reduced density matrix. 

Now the stochastic Schrodinger equation coupled to a finite bath, Eq. ([6]), 
will be implemented to calculate the energy decay, mean position and purity 
for the system particle inside two types of potentials, widely used in physics: 
the harmonic potential and the Morse potential. We also consider the time 
evolution of the quantum average position. 



3 Harmonic potential coupled to oscillators 

In this section we consider a particle inside a harmonic potential coupled to 
a discrete bath of HOs. This problem has exact solutions |2] which can 
be used to compare some quantities calculated with the approach used here. 
The frequencies ux are distributed in the interval (1.1,2.1) and are picked 
randomly from a frequencies generator, but once picked they will be kept 
fixed during the simulations, even for different bath realizations. In this ap- 
proach, due to the finite number of HOs in the bath, the spectral density is 
determined numerically and is always structured for low values of (see [TH] 
for more details in the classical case). The properties of the bath depend on 
the frequency's distribution of the N HOs [18]. In the limit — oo we as- 
sume to have three distinct continuous distributions: ohmic bath (quadratic), 
sub-ohmic and super-ohmic. Each oscillator will have a distinct {x\,yx) pair, 
randomly chosen numbers with zero mean and deviation one. For each coeffi- 
cient Cn of Eq. on]) we consider a rectangular initial conditions: c„(0) = 1 + ■ z, 
with n = 1, . . . , Umax, where rimax is the last energy level considered. The first 
i^max = 15 energy's levels from the harmonic potential are considered, with 
energy En = hjj{n — |). Using u = k = 1.0 the ground state energy level 
equals 0.5. Thus our initial state is given by 



P = Tr, p^it) =Y.M[cnC*jM[cnK 




nm 
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1 TT'max — -LO 

1^0) = ^= E \<Pn), (12) 
y'^max n=l 

which is an entanglement of all 15 system levels. In order to obtain the energy 
time decay and the average position as a function of the number of oscillators 
in the bath, Eq. (|6]) is integrated by a fourth-order Runge-Kutta integrator 
[38] and the values of c„(t) are introduced in Eqs. (j9]) and (fTOj) . 

In the first part of this section the frequencies of the harmonic oscillators from 
the bath will follow a quadratic distribution, making the spectral density to 
be linear J{u! < Ucut) oc u when — t- 00. 500 initial conditions for the pairs 
{x\,yx) will be considered for all the studied cases. To render the comparison 
easier the coupling strength will be kept constant, namely g\ = 0.01. 




Fig. 1. The average energy of a system, Eq. ([9]), coupled to a bath of one HO, with 
the coupling strength g = 0.01. The two insets show the corresponding average 
positions, Eq. ljlOp . for larger (left) and smaller (right) times. For comparison the 
circles (right) show exact results obtained from [2]. 

In Figured] the time evolution of the average energy (E), Eq. ([9]), is plotted 
for a bath containing one harmonic oscillator with the frequency u = 2.09. As 
inset graph is displayed the average position (q) from Eq. (fTO]) . One can clearly 
see that for times td ~ 170 the system loses around 20% of its initial energy 
for the bath, but is able to regain it at a later time t^j ~ 320. This exchange 
of energy between system and environment will repeat itself, with frequency 
around u ~ 0.037, for all integrated times (we checked it until t = 5000). 
Thus, the mean energy over the return times (where the energy transferred 
into the bath returns to the system) is almost constant. The return times of 
the energy are strongly dependent on the frequency of the bath oscillator but 
independent of the coupling intensity. This was checked for many coupling 
intensities. The time behavior of the energy was also confirmed by the time 
evolution of the average position: the system symmetrically oscillates very 
fast around the center (q) = 0, with an amplitude which follows the pattern of 
energy decay. Namely, the minimum of the energy corresponds to the minimum 
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Fig. 2. (Color online) The time evolution of all 15 energy levels £'„|c„(t)p 
of a system coupled to a bath of one HO. 

of the highest amplitude. Figure [2] shows the time evolution of all the 15 
individuals energy levels i?„|c„(t)p of the harmonic potential for the same 
bath used in Fig. [1] Numbers denote the energy quantum number n. One can 
clearly notice that while the lowest states remain practically unchanged, the 
higher energy levels oscillate following the energy exchange observed in Fig. [H 
Thus only the highest states are responsible for the energy exchange. 
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Fig. 3. (Color online) The average energy decay, Eq. ([9]), for a bath of 10 HOs, 
with the coupling strength equals to 0.01. By symbols we display the best fits. The 
two insets show the corresponding average positions, Eq. (fTO]l . for larger (right) and 
smaller (left) times. For comparison the circles in the left inset show exact results 
obtained from [2]. 

In Figure |3] the energy decay is shown in a semi-log plot for a bath of = 10 
HOs. In the time limit studied here (tmax = 500) one can clearly see that 
the system energy will not be regained so that the time average from the 
system energy is not constant anymore. However, we expect that for some 
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later times t > t^ax the energy will return to the system, since it is a fi- 
nite system. The qualitative behavior of the energy decay rate changes for 
different time intervals. We show by symbols the best fits of the energy de- 
cay, namely an exponential fit, a = 0.018, for short times where a Markovian 
dynamics is expected and a power-law fit, /3 = 1.245, for large times, where 
a non-Markovian dynamics occurs. We say that the qualitative change from 
exponential to power-law decay occurs close to times ta-^i3- As inset graph we 
display the corresponding average position. 




t 

Fig. 4. The time evolution of all 15 energy levels £'n|c„(t)p of a system coupled to 
a bath containing = 10 HOs. 

The time evolution of all the 15 energy levels of a system coupled to a bath 
which has = 10 HOs is displayed in the semi-log plot of Fig. HI Again it is 
clear that the higher energy levels decay faster. Something similar was shown 
to occur with the decoherence rate which is correlated to the highest occupied 
states [33] • In the long times region the energy levels are totally inverted: the 
highest energy level becomes the smallest one while the ground state energy 
becomes the most predominant one, followed by the second energy level and 
so on. By comparison with Fig. [3] we observe that the energy levels start to 
be totally inverted for times close to 250, i. e. ta^p where the transition from 
exponential to power-law decay occurs. At these times it can be observed in 
Fig-Hlthat the decay ratio dEn/dt decreases, meaning that the energy inserted 
in the bath starts to influence back the system dynamics, generating the non- 
Markovian behavior and the observed power law decay. 

The energy decay and the average position for a bath of = 20 harmonic 
oscillators is displayed in Fig. The system continues to behave dissipatively 
for the integrated time, with an exponential decay, {E) oc exp(— 0.030t) for 
< t < 50 and a power-law decay {E) oc t~^-^'^'^ for much larger times t > 230. 
For an intermediate time region, starting at t ~ 75 and ending at t ~ 190, an 
additional exponential decay, {E) oc exp(— 0.0226t), was observed. Such kind 
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Fig. 5. The average energy and its best fits for a batli of = 20 harmonic oscillators. 
As inset we shown the average position. 

of behaviors with two exponential decays was encountered until ~ 50 and 
completely disappeared for > 100. By keeping the same line of arguments 
of Fig. m we can say that, for a bath of 20 < < 50, the total inversion of 
the energy levels in times induces a power-law decay and the non-Markovian 
behavior. 



Fig. 6. The energy decay for a bath with N = 1,10,20,50,100,500 HOs. For 
N = 10, 20, 50 we observe exponential and power-law decays. 

In Figure|6]we summarize results of the energy decay for = 1, 10, 20, 50, 100, 
500 HOs. The first observation is that the energy decays faster for larger A^. 
For values A^ = 10, 20, 50 we observe an exponential decay for shorter times 
and power-law decay for larger times. Details of the decay exponents will be 
discussed later in Figs.[9]and[T0l As A^ increases the time interval of power-law 
decay decreases and vanishes totally for A^ = 100. Therefore, a non-Markovian 
motion is expected only for intermediate values of A^. For A^ > 50 a fiuctuating 
behavior is observed and will be explained below. 
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Instead of looking at the dissipation, it is also possible to analyze the purity 
from Eq. ( ITT]) which gives informations about the decoherence occurred in 
the system. For our initial entangled state from Eq. fll2p . the initial reduced 
density matrix PsiO) has 15 diagonal terms and 210 off-diagonal terms. The off- 
diagonal terms are characteristic of the quantum coherence. Total decoherence 
of the initial state occurs at times tdec when all off-diagonal terms vanish. For 
these times a statistical mixture of the 15 diagonal terms is reached and the 
the purity of the reduced density matrix is P{tdec) = 1/15 ~ 0.067. In this way 
it is possible to compare total decoherence times tdec with dissipation times. 




t 



Fig. 7. The purity for a bath with iV = 1, 10, 20, 50, 100, 500 HOs. 

In Fig. [7] we plot the time evolution of the purity, Eq. ( ITT]) , for the same 
parameters shown in Fig. O For = (no bath) the purity keeps constant 
to 1. For = 1 the purity diminishes and increases a small amount in time, 
but no total decoherence is observed. For = 10 we have one exponential 
decay with decay exponent ~ 0.024 for times t <67 and another exponential 
decay with exponent ~ 0.019 for times 75 < t < 190. For later time the 
purity decay obeys a power-law decay, as for the energy decay. In this case the 
total decoherence is observed for tdec ~ 140 when the purity crosses the long- 
dashed straight line marking the value P = 0.067. For the = 20 situation 
one can again notice three decay regimes: two exponential decays with different 
exponents 0.047, 0.018 (but they occur for the same time intervals as in A^ = 
10) and a power-law decay. The total decoherence is reached for tdec ~ 67. 
For the purity, the fluctuating behavior for long times starts with A^ = 50. By 
increasing A^ the fluctuation's amplitude decreases. 

Now, we explain the wavy pattern noticed for long times region t > 100. This 
phenomenon becomes more visible in the semi-log plot for higher values of A^ 
and it is a combination of two facts which act together. The first one is that the 
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difference between the ground state energy level and the second energy level 
increases by increasing A^, thus the ground state level becomes more and more 
predominant, compare Figures [2] and HI The second one is that the contribution 
of the ground state level to the energy exchange with the bath increases by 
increasing A^. For low values of the ground state energy is almost constant, 
while for > 50 the ground state level gives part of the energy to the bath 
and gain it back at a later time, creating this fluctuating behavior. The minima 
of such fluctuations, occurs close to times t = 60, 190, 310, . . ., which are very 
close to the times where a qualitative change in the energy and purity decay 
is observed in Figs. M and [71 



Figure [H] shows the phase space dynamics for one realization of the bath and for 
different values of A^ = 0, 1, 10, 20, 50, 100, 500, 1000. For A^ = 1 the evolution 
occurs over an ellipsis with a small width. When A^ increases to 10, 20 this 
width increases. It appears because the particle continuously loses some energy 
to the bath, but regain it later. See Fig. M^c) for A^ = 10 the particle motion 
just for initial times. At t = the particle starts on the outer side of the 
ellipsis. As times goes on, part of its energy is transferred to the bath and it 
moves towards the inner side of the ellipsis. Since there is a continuous energy 
exchange between system and bath, the particles moves continuously between 
the inner and outer side of the ellipsis. As A^ increases to 100, 500 the strong 
dissipative character of the system appears clearly. Particle moves towards the 
origin in phase phase, and stays there for all integrated times. 
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Fig. 8. Phase space dynamics for A^ 
right and top to bottom). 



0,1,10,20,50,100,500,1000 (from left to 



Now we turn our attention to other types of distributions, we consider the cases 
for which the spectral density has the form J{oj) = u'^, with < s < 2 [1]. 
If s = 1 the quadratic frequency distribution of the ohmic bath is recovered. 
For s < 1 we have the sub-ohmic bath, where lower frequencies in the chosen 
domain have a bigger contribution when A^ is sufficiently large. For s > 1 we 
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have the super-ohmic case, where higher frequencies appear more frequently. 




01^ , , — , — 1 < — < — — I 

10 100 1000 
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Fig. 9. Exponential fits, a, of the energy decay in the low times region for ohmic, 
subohmic, and superohmic baths function of A^. 

In Fig.[9]we plot the a exponents of the short time exponential decay for all the 
three types of baths: ohmic {s = 1.0), sub-ohmic [s = 0.1), and super-ohmic 
(s = 1.9). One can notice very small differences from one type of the bath to 
another. For a better visualization we plot the results in a log- log plot. One 
can clearly notice two distinct behaviors, shown in the figure by continuous 
lines. For low values of the exponent follows a power-law with the exponent 
equals to 0.749 while for high values of the best fit is another power-law 
function with the exponent equals to 0.362. 



o s=1.0 
□ s=0.1 



10 100 

N 

Fig. 10. Power-law fits, /3, of the energy decay in the long times region for ohmic, 
subohmic, and superohmic baths as a function of N. 

In Figure [To] we plot the /3 power-law exponents of the long times energy decay 
for ohmic, sub-ohmic, and super-ohmic baths. As before, there were noticed 
small differences in the exponents of different bath's type. We notice a fast 
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decrease of /3 when increases. For large A^, i.e. > 500, it is difficult to 
say that the energy decay follows a power-law function, due to the strong 
fluctuating behavior. 



4 Morse potential 



The Morse potential [39|40f41] is given by 



\/(r) = Z^Jl -e-"('"'*=f , (13) 



where r is the distance between atoms, r^, is the equilibrium bound distance. 




r 



Fig. 11. Morse potential for = 30, a = 0.08, rg = 0. With horizontal lines are 
depicted the energy levels in the r interval [—7.4,20], with continuous lines are 
shown the middle energies levels used in our simulations (see the text for details). 

The energy levels and the position elements rnm were determined by using 
the numerov method [42]. The results for the energy levels were comparable 
with the analytical form [39] : 



. l^e, I, [W^(" + i)]^ n4^ 

^«V^("'+o) —TR • (14) 



m ' 2' AD, 



e 



Figure [TT] shows in thick continuous line the Morse potential, given by Eq. 
(fT3l) . with the following values of the constants [40]: = 30, a = 0.08, = 0. 
With horizontal lines are depicted all the 38 energy levels, as calculated from 
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the numerov method. For the simulations we used as initial condition for the 
coefficients c„ a Gaussian wave packet form 0(r) = Cexp — '■'^"2^'' , where 
(7 = 3 and C is the normalization constant. Thus, we will have a wave packet 
with the peak at the energy level 16, placed in the middle of the energy 
spectrum, corresponding to the values 5.03 to 12.32, depicted by continuous 
horizontal lines in Figure [TTJ Here, the coefficients c„ have the real Gaussian 
form defined above and they fulfill the equation I]n=9 |cnP(0) = 1. 




Fig. 12. The energy decay, Eq. ([9]), for Morse potential. The number of HOs in the 
ohmic bath is varied: 0, 1, 10, 50, 100. 

The energy decay, Eq.Q, for a bath having iV = 0, 1, 10,50, and 100 HOs 
is plotted in Fig. [121 It was considered that all of them have the same mass 
m = 1 and experience the same coupling strength with the system, gx = 
g = 0.001, which is one order of magnitude smaller when compared to the 
HO case considered before. Here, the frequencies of the oscillators follow a 
quadratic distribution in the interval [1.1, 2.1], and for each bath, once chosen 
the frequencies, are kept constant. For each case, namely a bath having = 
1,10,50, and 100 harmonic oscillators, we will vary only xx and yx, which 
are real distributed numbers with zero mean and deviation one. For all the 
non-zero cases {N 7^ 0) we took 500 different initial conditions {xx,y\), except 
the case of = 100, where due to the high CPU time we average over 200 
initial conditions. For the simulations we consider all the energy levels of the 
system, described by a wave packet of 38 energy levels, e„, corresponding to 
38 wavef unctions, 0„ (and implicitly the c„ coefficients). 

It can be easily seen that in the case of no bath = the total energy is 
conserved in time. When we have the case = 1, in the chosen time interval, 
the total energy decays slowly but returns to the system at t ~ 150. For 
A^ > 10 HOs the energy decays exponential in time, with some fluctuations, 
due to the stochastic average. The power-law decay observed in the HO case 
for later times is not seen here since the coupling to the bath is much smaller 
and decays occur much slower. In fact the purity (not shown) also decays much 
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slower when compared to the HO case. The exponential energy decays can we 
0.01^ ' — ' — 1 




Fig. 13. The exponent a of the energy decay as a function of N. 

summarized in Fig. [131 where the values of the exponential coefficient a of 
the energy decay are plotted as a function of the number of the oscillators of 
the bath, A^. For a better visualization the results are shown in a log- log plot. 
The situation is different than the power-law fits encountered for the harmonic 
potential. Here, we found that the best fit is a quadratic function given by 
10-^[0.0307A^2 ^ 2.25iV + 22.9]. Thus exponential decays increase with A^. 



5 Conclusions 

In this paper we use the stochastic Schrodinger equation for zero temperature 
to study the energy decay for a system particle, situated in a harmonic poten- 
tial and in a Morse potential. The system is coupled to a bath composed of 
a finite number N of uncoupled harmonic oscillators. For the discussions we 
chose, in the limit N — )■ oo, an ohmic bath distribution, but we also study the 
cases of sub-ohmic and super-ohmic baths. 

In the case of a harmonic potential, with coupling intensity to the bath 7 = 
0.01, it was observed that for very small numbers N the energy is exchanged 
back and forth into the bath. For intermediate values of N around 10 < < 
20, the time average energy starts to decay, transferring partially its energy 
to the bath and we do not see any energy regain for the integrated times. For 
these values of N we observed two exponential energy decays (with different 
exponents) for small times and a power-law decay for large times. For relatively 
higher values of (> 50) a strong exponential decay was observed and no 
power-laws. Therefore, non-Markovian dynamics is expected for intermediate 
values of A^ (10 ^ N < 20) and a Markovian motion for larger values of A^. 



15 



The exponents for the exponential decays increase with obeying two power- 
laws. We analyzed also the purity of the system which gives informations 
about the decoherence process. Essentially decoherence occurred in all cases 
for short times, where exponential decays were observed for both, the energy 
and purity. For a system situated in a Morse potential, with coupling intensity 
to the bath 7 = 0.001, the high computational effort forced us to restrict our 
simulations to shorter time evolutions. Due to such small couplings, the energy 
and purity decays occur much slower when compared to the harmonic case and 
we just observed the exponential decays. Different from the harmonic case, the 
exponents of the exponential decay increase quadratically with A^. 

One of the main goals of this work is to show general features which occurs in 
systems coupled to an increasing number of bath constituents. Therefore, in 
addition to the above main general results, we would like to conclude with some 
relevant characteristics observed in the simulations but not mentioned along 
the text: The energy decay rate increases with increasing coupling strength 
between system and bath, but the time for energy regain is independent of the 
coupling. The energy regain time depends on the mean bath frequency and on 
A^. For larger results are totally independent on the numerical frequency 
generator. For smaller values 1 < A^ < 50, the values of exponential and 
power-law exponents for the energy decay may change for different frequencies 
generator, but the overall qualitative dynamics is equivalent. 
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